The oxygen isotope effect on critical temperature in superconducting copper oxides 
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The isotope effect provided a crucial key to the development of the BCS (Bardeen-Cooper-Schrieffer) 
microscopic theory of superconductivity for conventional superconductors. In superconducting 
cooper oxides (cuprates) showing an unconventional type of superconductivity, the oxygen isotope 
effect is very peculiar: the exponential coefficient strongly depends on doping level. No consensus 
has been reached so far on the origin of the isotope effect in the cuprates. Here we show that the 
oxygen isotope effect in cuprates is in agreement with the bisoliton theory of superconductivity. 

PACS numbers: 74.62.-c: 74.62. Dh; 74.72.-h 



Just before the development of the BCS theory 
(Bardeen-Cooper-Schrieffer) [1], it was experimentally 
found that different isotopes of mercury have different 
critical temperatures, and TcM 1 / 2 = constant, where M 
is the isotope mass [2,3]. In the framework of the BCS 
theory, the expression for the critical temperature T c is 
given by fcsT c = 1.14 Tuvd exp(— 1/A), where u>n is the 
phonon frequency at the edge of the Debye sphere; A is 
the electron-phonon coupling constant (= N(0)V); ks 
is the Boltzmann constant, and h is the Planck constant 
divided by 2ir. So, the BCS theory has provided an expla- 
nation for the isotope effect because the Debye frequency 
varies as a Af -1 / 2 , and consequently, the product 
T C M X I 2 is constant for a given superconductor. For the 
majority of superconducting elements, the exponential 
coefficient a is indeed close to the classical value of 1/2. 
Later on, it was discovered that the situation is more 
complicated than it had appeared to be. For some con- 
ventional superconductors, the exponent of M is not 1/2, 
but near zero. For example in Ru and Zr, a = 0±0.05. 
Thus, the isotope effect is not a universal phenomenon, 
and can be absent even in conventional superconductors. 
Nevertheless, this effect played the decisive role in show- 
ing the way to the correct theory of superconductivity in 
metals. 

In superconducting cooper oxides (cuprates) [4], the 
situation is very peculiar. The oxygen isotope effect was 
found to be extremely small in optimally doped cuprates. 
This fact was initially taken as evidence against the BCS 
mechanism of high- T c superconductivity (true) and, mis- 
takenly, against the phonon-pairing mechanism (false). If 
the pairing mechanism is different from the BCS mecha- 
nism, this does not mean that phonons are irrelevant. 

In fact, there is a huge isotope effect in the cuprates. 
Figure 1 shows the oxygen ( 16 vs 18 0) isotope-effect co- 
efficient olq — d(lnT c )/d(lnM), where M is the isotope 
mass, as a function of doping level in La2-zSr. E Cu04 
(LSCO), YBaaCugOg+x (YBCO) and Bi 2 Sr 2 CaCu20 8 
(Bi2212). In the plot, one can see that the oxygen-isotope 
effect in the cuprates is not universal: it is system- and 
doping-dependent. In the undcrdoped region, ao can 
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FIG. 1. Oxygen isotope effect ( 16 vs 18 0) in LSCO, 
YBCO and B12212 (taken from Ref. [5]). 



be much larger than the BCS value of 0.5 (according 
to the BCS theory, the isotope effect cannot be larger 
than 0.5). In the optimally doped region (p ~ 0.16), the 
oxygen-isotope effect is indeed small. With the excep- 
tion of one point in LSCO, p = 0.125 (the so-called | 
anomaly in LSCO), the doping dependence of the coef- 
ficient ao is universal in these cuprates. The coefficient 
ao has a maximum at p — > 0.05. 

Even if a few attempts have been made in the literature 
to explain such a peculiar isotope effect in the cuprates 
[6], no consensus has been reached so far. The main pur- 
pose of this paper is to show that the isotope effect in 
cuprates can be well understood in the framework of the 
bisoliton theory of superconductivity for cuprates [7-11]. 
One of the chief results of the bisoliton model is that the 
potential energy of a static bisoliton, formed due to local 
deformation of the lattice, does not depend on the mass 
of the elementary lattice cell. This mass appears only in 
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FIG. 2. Critical temperature as a lunction of Debye tem- 
perature for different cuprates [12]. 



the kinetic energy of bisolitons. It is a paradox: the ef- 
fect of isotope substitution on the transition temperature 
manifests itself when the electron-phonon interaction is 
weak, and can disappear when the electron-phonon in- 
teraction becomes stronger! 

The cuprates are materials with strongly correlated 
electrons. The Fermi-liquid approach is not applicable to 
this class of materials because the position and motion 
of each electron in these materials are correlated with 
those of all the others. If in metals, the electron-phonon 
coupling is week and linear, in materials with strongly 
correlated electrons, it is moderately strong and nonlin- 
ear [11]. Therefore, phonons interact with charge carriers 
in conventional superconductors and in the cuprates in 
a different manner. For example, in BCS-type supercon- 
ductors, T c increases with lattice softening, while in the 
cuprates, T c increases with lattice stiffening, as shown 
in Fig. 2. In the BCS model of superconductivity, only 
acoustic phonons participate in electron pairing, while 
tunnelling studies of the spectral function a 2 F(ui) show 
that in cuprates, the charge carriers are coupled not only 
to acoustic phonons but also to optical ones [11]. The 
function a 2 F(uj) is the parameter of the electron-phonon 
interaction in the Eliashberg equations, which character- 
izes the coupling strength between charge carriers and 
phonon vibrations. The energy of these optical phonons 
is about 73 meV. In the cuprates, the 73 meV branch is 
associated with half-breathing-like oxygen phonon modes 
that propagate in the Cu02 plane. 

The bisoliton model of superconductivity developed by 
Davydov for organic and high-T c superconductors [7-9] 



is based on the presence of one dimensionality in the sys- 
tem. So, in the bisoliton theory, the electron transport is 
considered along one-dimensional chain (when the bisoli- 
ton model was created there was no knowledge about 
charge stripes in the cuprates). Excess quasiparticles 
present on the chain can be paired in a singlet state due 
to the interaction with local chain deformation created 
by them. The potential well formed by a short-range 
deformation interaction of one quasiparticle attracts an- 
other quasiparticle which, in turn, deepens the well. The 
energy of a pair of quasiparticles moving in the field of 
local deformation of the lattice with a velocity v, called 
a bisoliton, is given by 



E(v) = E p (0) 



m bs v* 



(1) 



where E p (0) characterizes the position of the energy level 
of a static pair of quasiparticles beneath their Fermi level 
Ep , and is the effective mass of the bisoliton. The 
potential energy of a static bisoliton E p (0) is independent 
of the mass of the elementary cell M. At a small density 
of doped charge carriers, the effective mass of a static 
bisoliton can be represented as 



m bs 



8MA P 

2m + H^- 



(2) 



where m is the effective mass of quasiparticles; A p is 
the pairing energy gap; k is the coefficient of longitudi- 
nal elasticity, and a is the distance between atoms on 
the chain. The mass of a bisoliton exceeds two effective 
masses of quasiparticles forming a bisoliton. It is worth 
noting that the weakest point of the bisoliton model of 
superconductivity is that the Coulomb repulsion of elec- 
trons is not taking into account. In spite of this fact, 
the bisoliton theory gives pretty accurate values of the 
Cooper-pair size in hole- and electron-doped cuprates 
[11]- 

Substituting the expression for mi, s into Eq.(l), one 
obtains 



E(v) = E p (0) + mv 2 



AMA p v 2 



(3) 



In this expression, only the last term contains the mass 
of the elementary lattice cell M. Therefore, in the frame- 
work of the bisoliton model, the isotope effect is deter- 

A v 2 

mined by the product J 2 . Let us analyze its doping 
dependence. The divisor ha 2 is a weak function of p. 
The doping dependence of the pairing energy gap A p in 
hole-doped cuprates is shown in Fig. 3. Figure 4 depicts 
the doping dependence of the velocity of nodal quasipar- 
ticles, determined by angle-resolved photoemission mea- 
surements. In Fig. 4(b), one can observe that the nodal 
high-energy quasiparticles in hole-doped cuprates move 
quicker in the underdoped region than in the overdoped 
region. Assuming that nodal quasiparticles with high 
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FIG. 3. The phase diagram of hole-doped cuprates [13,11]. 
A c is the phase-coherence energy gap, and A p is the pairing 
energy gap. 



energy play an important role in the occurrence of su- 
perconductivity in the cuprates, we obtain that, in the 
framework of the bisoliton model, the product A p vf 1E 
which determines the isotope effect has a doping depen- 
dence similar to that of ao shown in Fig. 1. Thus, in 
the framework of the bisoliton theory, the doping depen- 
dence of ao in Fig. 1 is natural. We are not going to 
discuss here the | point in LSCO, as well as the weak 
uprise of ao (p) in the overdoped region which can be due 
to several effects. 

To conclude, we showed here that the isotope effect 
in the cuprates can be qualitatively understood in the 
framework of the bisoliton theory of superconductivity 
[7-11,5]. It follows from the fact that the potential en- 
ergy of a static bisoliton, formed due to local deformation 
of the lattice, does not depend on the mass of the elemen- 
tary lattice cell. In the framework of the bisoliton model, 
this mass appears only in the kinetic energy of bisolitons. 




FIG. 4. Velocities of nodal quasiparticles as a function of 
doping level p in hole-doped cuprates [14] . (a) The low-energy 
(< 50 meV) Fermi velocity, and (b) the high-energy (> 50 
meV) velocity. 
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